DNA as a one-dimensional chiral material: Application to the structural transition 

between B form and Z form 
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A dynamical model is presented for chiral change in DNA molecules. The model is an extension 
of the conventional elastic model which incorporates the structure of base pairs and uses a spinor 
representation for the DNA configuration together with a gauge principle. Motivated by a recent 
experiment reporting chiral transitions between right-handed B-DNA and left-handed Z-DNA [M. 
Lee, et. al., Proc. Natl. Acad. Sci. (USA) 107, 4985 (2010)], we analyze the free energy for the 
particular case of linear DNA with an externally applied torque. The model shows that there exists, 
at low temperature, a rapid structural change depending on the torque exerted on the DNA, which 
causes switching in B and Z domain sizes. This can explain the frequent switches of DNA extension 
observed in experiments. 

PACS numbers: 87.14.gk, 82.37.Rs, 87.15.A- 
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I. INTRODUCTION 

Recent advances in experimental techniques have made 
it possible to perform experiments on the mechanical re- 
sponse of DNA molecules involvingstretching and twist- 
ing of single DNA molecules [il,l2|. It has been ob- 
served that external forces can bring about changes in the 
macroscopic conformation as well as the base-pair struc- 
ture of DNA. For example, frequent structural changes 
have been observed between right-handed B-DNA and 
left-handed Z-DNA when minute negative torque is ex- 
erted on a DNA molecule Q. In the regime of larger 
force, other structures, such as S-DNA, P-DNA, etc., 
are found [4], and the coexistence of these structures in 
stretch- twist diagrams has also been reported Q. 

Along with the experimental studies, various theo- 
retical attempts have been made to describe the me- 
chanical responses of DNA. One typical theory known 
as the elastic rod model enables a description in terms 
of a simple mechanism of supercoiling and its effects 
[H, Other significant theoretical models, such as 

the Poland-Scheraga model [9|] and the Peyrard-Bishop- 
Dauxois model [10], are concerned with the problem of 
denaturation, and have been extended so as to include 
the coupling between DNA twisting and denaturation 
pj| . The interaction between DNA supercoiling and 
denaturation has been studied by extending the classi- 
cal elastic rod model [12j. An effective potential model, 
which predicts twist-stretch coupling in the mechanical 
response of B-DNA, was presented in [l^. Statistical 
models that phenomenologically describe various struc- 
tural transitions were developed in [l^. Several meso- 
scopic models have also been developed, which can de- 
scribe the interaction between DNA conformation and 
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FIG. 1: (Color online) (a) The DNA molecule consists of two 
sugar-phosphate chains bridged by base pairs, forming the 
double helix. Moving frame vectors (6, n, t) are defined along 
the helical axes parameterized by a continuous parameter s. 
(b) The frame vectors are specified by the depicted Euler 
angles (6>,0,x)- 



melting structural transition [TBI, |l6|. However, there has 
not been a mechanical model of DNA that describes the 
interplay between the global configuration and the vari- 
ous intrinsic structures of DNA base-pairs, including B, 
Z, S structures. 

The purpose of this article is to propose a phenomeno- 
logical model that combines the aspects of conformation 
and intrinsic structure to obtain a clear picture of DNA 
mechanics. Our model adopts the Landau theory, which 
consists of the free energies for helical structures of base 
pairs connected by hydrogen-bonds, elastic deformations 
of DNA configuration, and the interaction between them. 
The interaction is introduced by using the gauge princi- 
ple, with the spinor representation of DNA configura- 
tion. Monte Carlo simulation of the model is used to 
examine B-Z structural transitions of linear DNA un- 
der the application of external torques. It is found that 
there exists a cooperative effect, depending on both the 
temperature and the torque exerted on the DNA, which 
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can be described in terms of the probability to create 
kink-antikink pairs. This cooperative effect causes sharp 
switching in B and Z domain sizes, which can explain 
the frequent switches of DNA extension observed in re- 
cent experiments [3|]- 



II. MODEL 

We use two kinds of order parameters: one is the or- 
der parameter for the conformations of DNA and the 
other is that for the intrinsic base-pair structures, which 
are determined by the hydrogen bonds between them. 
We first construct a moving frame that describes the 
conformation of the double-stranded DNA of N base- 
pairs (bp's). Figure [TJa) shows a continuous parameter 
s representing the number of bases from one end of the 
molecule, where the DNA helical axis is parametrized as 
r{s). The normalized tangent vector along the axis is 
defined as t = ^ ^ , where v is the increase per bp along 
axis ||(ir/(is||. A moving frame is defined at each point s, 
by the set of three mutually orthogonal vectors (6, n,t), 
where b is the normalized vector in the direction of base 
position and n is specified as n = t x 6 to form a right 
handed coordinate system. 

We now introduce an order parameter to describe 
the configuration of the moving frame (6, n,t), which 
is represented by a spinor ^: ^ = ('0i,V^2) ^ with 
IV^iP + l'^2p = 1- The spinor is parameterized by a set 
of Euler angles, as in [17]. Using Euler angles (6>, 0, x) as 
depicted in Fig. [TJb), the spinor is written as 



^ = (cos-e 2 ,sm-e 2 y 



(1) 



and t can be written as an average of the Pauli vector a 
by the spinor: t = "^^a"^ = (sin 6> cos ^, sin6> sin^, cos 6>), 
where a is given by (a^, <j^, a^) with the following Pauli 
matrices: 
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^ = '0 -1 



We next consider an explicit form for the free energy of 
the DNA conformation, in terms of a functional of ^. To 
this end, we set the following requirements (i) rotational 
symmetry and (ii) inclusion of up to first derivative with 
respect to s. Then, the simplest form satisfying these 



criteria is given by F = Tds^ where 



ki d^t k2 



ds 



(2) 



With the use of the Euler angles, J^g reads f [(^)^ + 
{^fsm'e] + q{^cose + ^)^ with B = fci/4, C = 
(^1 + A:2)/4, which is the well-known free-energy density 
of an elastic rod with isotropic bending elasticity B and 
torsional elasticity C [l8|. Note here that, if an external 
force / in the z-direction and torque r are applied to 



the DNA ends, the free energy density becomes J^s — 
fvcos(0) — rLk^ where is the linking number of the 
DNA 0. 

Now we introduce the order parameter p represent- 
ing the internal structures of base pairs, via gauge cou- 
pling with the conformational spinor. Consider the gauge 
transformation ^ ^exp(ia). If a is constant, the in- 
variance of the free energy under the transformation is 
apparent, but if a has 5-dependence, the local gauge in- 
variance does not hold. To keep the gauge invariance [l9|, 
we define the order parameter p for the internal structure 
of base pairs as a gauge field, which we call chiral field, 
p is the rotational angle between successive base pairs 
in the helix, p > for right handed rotation and p < 
for left handed rotation. After progressing ds, we have 
X ^ X + pds and 



^11(5 + ^5) = {t/jie' 



-i^ds 



i^^ds, 



which just gives the parallel transport. Since the dis- 
placements of sugar-phosphate strings from the paral- 
lel transport ^"(s) gives rise to an elastic potential, we 
make the replacement oid/ds D/Ds in Eq. ([2]), where 
D/Ds is the covariant derivative such that 



lim 

ds^O 



"^{S^ds) -^ll(5 + c/5) _ 

ds ds 
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The resulting free-energy density becomes Ts + 

where Tsp = -^Cp(^t^ - + thereby we 

have a new coupling term Tsp between the spinor field 
and the chiral field. In terms of the Euler angles, the 
coupling term is given by 



-Cp{ 



dx^d^ 
ds ds 



cos 0) 



Cp^ 
2 



from which we see that is conjugate to the angular 
velocity with respect to s around t (1^3 = ^ + ^ cosO). 
Combining the above two terms, the total free-energy 
density J^s+p can be written in terms of the Euler angles: 



c , 



(3) 

This is similar to the elastic rod model [6], but differs in 
the point that there is a couphng with the chiral field p. 
The conventional elastic rod model does not allow the in- 
clusion of such a coupling term, because it does not take 
into account the deformation degrees of freedom coming 
from the base pairs, whereas in our case the spinor repre- 
sentation enables us to take into account such degrees of 
freedom naturally as a gauge field through the covariant 
derivative. 

Now Landau's free energy for base-pair structures is 
defined by a general functional of the above-introduced 
chiral field p{s) and the length scale v{s) = \\dr/ds\\. The 
free-energy density is given by 



di.dp.2 



-(-) +-. + Vip,v,...). (4) 
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Here the last term is the effective potential, whose local 
minima, (pb,^b,...)^ (pz,^z,...)^ (Ps,^s,...)^ etc., can 
be assigned to the locations of the base-pair structures 
B-DNA, Z-DNA, S-DNA, etc., respectively. On the other 
hand, the gradient terms represent additional free-energy 
costs due to changing base-pair structures along s. 

In this way, we have constructed the total free-energy 
density for DNA as 

T = Ts^p^Tb' (5) 

It should be noted that the coefficients 5, C can be ex- 
tended to be functions of the order parameters for base- 
pair structure [20| . Our model contains the conventional 
elastic model, as a limiting case of the rigid base-pair 
limit, V{p) = —S{p — ujq). Namely, the latter is obtained 
by substituting ujq foi" P in Eq. (j3j). In the rigid sugar- 
phosphate limit C ^ oo, it recovers the model given in 
Ref. [21| , becoming nonlocal in the sense that it contains 
the second derivative of the order parameter. 

III. CHIRAL TRANSITION INDUCED BY 
EXTERNAL TORQUE 

As the first application of our model to realistic con- 
ditions, we consider the experiment performed by Lee, 
et.al., in Ref. [3]. In the experiment, shown schematically 
in Fig. [21 the ends of a DNA molecule were fixed with a 
glass cover slip and a magnetic bead, respectively. By ap- 
plying negative torque to the DNA with the bead, minute 
negative superhehcity was induced in the DNA and then 
an inter-winding DNA configuration, called plectoneme^ 
was formed. Below the plectoneme, in the linear part 
of the DNA, there was a poly(GC) sequence that tends 
to change structures between B-DNA and Z-DNA. They 
measured frequent sharp changes of the DNA extension 
between the cover slip and the bead, as well as frequent 
B-Z transitions in the poly(GC) part. Therein, the B-Z 
transitions induce winding and unwinding of the plec- 
tonemic conformation, resulting in the observed sharp 
switches between small and large DNA extensions, re- 
spectively. 

To analyze this experiment, it is sufficient to consider 
a model of a linear DNA molecule with B-Z structural 




FIG. 2: Schematic representation of the experiment in Ref. J3\ 
(see text). 



transitions. By setting 6> = 0, ^ = in Eq. (j3j), we obtain 
the total free-energy density given by 

where we adopt the following values: C = 91 x 10~^^J, 
di = Di/uol with Di = 4.1 x 10"^^ J and uq = 0.6 
rad/bp, and r is the torque applied to the end of the 
linear DNA molecule. The potential density V{p) is given 
by 

V{p) = Vo[{p/ouof-l]'^rcP, (7) 

with Tc = -7.9 X 10-21J [22j and = 6 X IQ-^Oj, 
where the B-DNA and Z-DNA conformations correspond 
to the minima at around p = ujq and p = — cjq, respec- 
tively. These parameter values are determined as fol- 
lows. We expect the twist persistence length is 75nm, 
which gives C 75nm x kBTn/OMnm 221kBTR 
91 X 10~^^J, where Tr denotes room temperature and 
ksTn = 1.38 x IQ-^^JK-^ x 300K = 4.1 x IQ-^^J 
is used. In our model, ^{dp/ds^ gives the energy 
cost for structure change, which was already given in 
Ref. [21] by 2kBTR = 8.2 x IQ-^^J. From the equal- 
ity ^(2(jJqY — one obtains d\ = Di/uj^^ where 
Di = ksTn = 4.1 x 10"^^J. As to Vb, we adopt the 
value Vo = 6x 10-2°J(- IbkBTn) according to Ref. [2l|. 
With this value, the domain wall energy of our model 
{E^_z = 8/3^/2diVoUJo [23|) agrees with the "B-Z junc- 
tion energy" 5 Kcal/mol, that was estimated in Ref. (2^ . 
The free energy of Z conformation is Ez = 2.3/cbTr 
thereby Tc satisfies 2ujoTc = — 2.3/cbTh, which gives 
Tc = —7.9 X 10~^^. Note that we have chosen the depths 
of minima that correspond to normal physiological condi- 
tions [25|. Furthermore, for simplicity, we have neglected 
the 'U-dependency supposed in the general expression of 
Eq. gj. 

To perform numerical simulation, we discretize the pa- 
rameter s. The resulting free energy is given by 

N 

F = Y.^V{p,)^-{x^-X^-l-p^? 

+Y(pi - pi-if - ^ixi - Xi-i)], 

with Xo = and po = ^o- By using the replica ex- 
change Monte Carlo method [26|, we sampled the equi- 
librium states, under various conditions of the external 
torque r and temperature T. With every step, this 
method updates all replicas Xr with different temper- 
atures Tr (r = 1,2,...) simultaneously and then ex- 
changes these rephcas, to enhance the ergodicity of these 
samples. In the simultaneous updates, each replica = 

{p^;\/^xf.P^2\^Xp.---JN^ where Axf^ ^ 

xt^ - xt-i (i-e-' xt^ = Y!j=i ^X^P) follows the follow- 
ing steps: 

Step 1: Draw a random integer i from 1 to N . 
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FIG. 3: (Color online) Exchange Monte Carlo results for N — 
100 with 1000 samples: {di)xn/ri (n — 50), (h)wz, (c)Log 
plot of Wz, (d) Log-Log plot of wz, as functions of r — Tc for 
/3 = 1(#, red), 3(B, yellow), 5(*, green), 10(A, blue), and 
14(T, black) in units of Vb(= 6 x 10~^°J), where the vertical 
error bars on data points represent their variances, wz is 
symmetric with respect to {r,wz) — (tc, 0.5). The i-ih bp is 
considered to be Z-structure if yo^ < 0. Solid lines in (d) are 
theoretical curves of Eq.([T5)) with fitted values Eb-z = O.llVo 
and 7 = 0.38 [13]. Between x(^) ^ind wz, the approximate 
relation dx/ds ^ ps + (pz — Pb)wz -\- r/C holds, where pB — 
0.6, pz = -0.6. 



(r) 



Step 2: Ax\ 

~(r) 

Pi = P. 

where e 



and p, 



(r) 



^^^+e, R 
and ev 



are supposed to be shifted to 
and Axf^ = Axf^+e^ 
are the values of maximum shifts 



for p and Ax, respectively, and R- 
number between —1 and 1. 



1,1 



is a random 



Step 3: Compute the free-energy change SF = F{Xr) — 
F{Xr) and W = exp(-/3^^F), where /3r = 




Wz 



(a) Low temperature 



T<TC ^<C—TC T>TC 




(b) High temperature 



FIG. 4: Schematic illustration of torque responses of wz. (a) 
At low temperature, wz behaves like switch. Three curves 
represent there cases of torque distributions: the average r 
satisfies t < Tc (left), r = Tc (middle), and r > Tc (right). The 
fluctuations 6t are given by Eq. jS]). Shaded regions in the 
distributions represent the ratios of Z-DNA conformations, 
(b) At high temperature, wz changes gradually as a function 
of T, compared to Sr. Hence, the ratios of Z conformations 
also change gradually, and do not show switching behaviors. 



(kBTr) ^ with Boltzmann's constant ks- 
Step 4: If W > i?o,i then make changes Ax-^^ = Ax-^^ 



and p. 



(r) 



(r) 

p\ ; otherwise do not change. 



After these simultaneous updates, for a randomly chosen 
neighboring pair of replicas, say Xj. and X^+i, exchange 
their configurations by the following acceptance criterion: 

Step 1: Compute the cost function 

A = {pr+l - Pr)[F{Xr) - F(X,+i)] 

and W = exp(— A). 

Step 2: If W > Ro,i then exchange replicas as 
{Xr^ Xr-\-i) = (X^+i,X^); otherwise do not ex- 
change these replicas. 

In the simulations, we set {Pr} = {0.01, 0.03, 0.05, 0.07, 
0.1, 0.3, 0.5, 0.7, 1, 1.25, 1.5, 1.75, 2, 2.5, 3, 3.5, 4, 4.5, 
5, 5.5, 6, 6.5, 7, 7.5, 8, 8.5, 9, 9.5, 10, 10.5, 11, 11.5, 12, 
12.5, 13, 13.5, 14}, and confirmed the frequent exchanges 
of replicas. All the states are sampled at every lOOA^ 
exchange Monte Carlo steps, after equilibration. 

Fig. [3l^a) shows the average twists as functions of r. 
We see that, at high temperatures compared to Vb, the 
response to r is approximately uniform. In contrast, at 
relatively low temperatures, xa^ are susceptible to small 
differences in r especially at around r = Tc. This result 
is in accord with the experiment of Ref. [[3] . Namely, in 
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the experiment, the torques apphed to the hnear part 
of DNA can not be controhed to be a certain value but 
actually fluctuate around the value as 



where 



Sr - ^CksT/N - 0.032Vb, 



(8) 



due to the disturbance from the rest of the DNA and/or 
the physiological buffer, and thus xa^ changes widely 
due to the high susceptibility to r at around Tc, since 
the experiment was carried out at low temperatures 
(/3Vb ~ 14). The large changes in xn alternate the wind- 
ing number of the plectoneme, resulting in the change of 
DNA extension observed in the experiment. 

We shall now confirm that the sensitivity of xn to r 
is caused by the B-Z transition. For this, the Z-DNA 
ratio, wz^ is plotted as functions of r in Fig. [3l^b)-(d). 
This clearly shows that the sharp sensitivity of xn to r 
at low temperatures found in Fig. [3fa) is coincident with 
the sharp decrease of wz^ while the mild dependency of 
Xn on r at high temperatures is accompanied with slow 
decrease of wz- Hence, our model clearly shows that the 
experimentally found sharp changes can be attributed 
to the dependency of B-Z transitions on torque, which 
becomes singular in the low temperature limit. 

Let us discuss the details of the structural transition 
induced by external torque and the resulting singular re- 
sponses, shown in Fig. [3l In low temperature such that 
T /Vq ~ (room temperature satisfies this condition), 
the response to torque is sharp, and then behaves like 
a switch. Hence, for example, if r — Tc > is suddenly 
changed to r — Tc < 0, then all DNA structures nucleate 
from false ground state (B-DNA) to true ground state (Z- 
DNA). This point is illustrated in Fig. HJ^a). Otherwise, 
if the superhelicity is set around "the midpoint value," 
then the torque fluctuates at around r ~ Tc, which in- 
duces the rapid structural switching observed in Fig. 5A 
in Ref. Namely, we have confirmed that this switch- 
ing induces the changes both in wz and in the change of 
DNA extension, dL. In contrast, when the temperature 
is high (T/Vb > 1), the response becomes gradual, as 
shown in Fig. IH^b). In this case, we can neither expect 
the switching behavior in wz nor dL. In this way, with 
the use of our model, we have succeeded in clarifying 
the condition for stepwise change and for interconversion 
between these states observed in Ref. ^3|] . 



IV. NUCLEATION THEORY 

Our findings can be explained in terms of the proba- 
bility of kink-antikink nucleation and the average size of 
the metastable domain p7j . 

The variational equation 5F = reads 



dx 
ds 



T 

c 







and 



ds^ 



dVM 
dp 



(9) 



(10) 



VM = Vo 



-\ 2 



^ -1 



(r - T,)p. (11) 



Vef[ is depicted in Fig. [5fa) for various r. For r > Tc, B 
conformation is stable, while Z conformation is unstable. 
In contrast to this, for r < Tc, Z conformation is stable, 
while B conformation is unstable. At r = Tc, KfT is 
symmetric and these conformations are equally probable. 
The probability P to form the nucleation is given by 



[1 



p/3£^nucl(r)l 



(12) 



where E^iudi^) the nucleation energy at torque r. Since 
the nucleation starts with creating kink-antikink pairs, 
Enuc\ is roughly estimated to be twice the domain-wall 
energy £;dom U^ '- 



(13) 



Then, the torque-dependence of the domain wall energy 
is assumed to be 



Edom{r) = Eb-z + 7k - 



(14) 



where Eb-z is the domain wall energy at r = Tc, shown in 
Fig.[5][b), and 7 is the coefficient of torque in the domain 
wall energy, shown in Fig. [5l^c). 

The distance between centers of neighboring 
metastable domains is estimated to be 1/P. Once 
a metastable domain is created, it grows to a size r with 
the probability exp(— /3AVr)/A/', where AV is the free- 
energy density difference and A/" = J^^^ exp(—/3A Vr)(ir. 
The average domain size is estimated as 



(r) 



i/p 



-(3AVr 



= P 



-\- X — expx 
x{l — expx) 



with X = PAV/P. From these results, the ratio of do- 
main size between the metastable and the stable domain 
is given by (r) : P~^ ~ (^)- If > ''"c, then B-DNA is 
stable and Z-DNA is metastable and thus wz = {r)P. 
Otherwise, if r < Tc, then Z-DNA is stable and B-DNA 
is metastable and wz = 1 — (r)P. Putting these estimates 
together, we obtain 



1 



Wz 



x{l — e^) 



(15) 



with 



^ /9(Vz - Vfe)/P, - Vb = 2wo(r - t,). (16) 



Let us confirm that Eq. (fT5]) describes our numerical re- 
sult of Fig. [3l The solid curves in Fig. [3jd) plot this re- 
lation (p!5]) . which agree with the numerical calculations. 
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FIG. 5: (Color online) (a) Illustration for Veff in Eq. (pT]) as functions of p\ for r < Tc (red), r = Tc (black), and r > Tc (blue). 
Veff is symmetric at r = Tc. Figs, (b) and (c) schematically show the domain- wall energy given by Eq. ([14]): (b) Eb-z is the 
domain wall energy for the symmetric r = Tc case. In this case, the domain-wall energy corresponds to the shaded region, (c) 
For T / Tc case, the domain wall energy requires the excess energy, shown in darker shade of gray. This figure shows that the 
excess energy is proportional to |r — rc|, as in Eq. ()14|) 



We can also explain the singularity. Expanding wz 
around r = Tc gives 

wz - 0.5 - P-^/3uJo{r - rJ/G, (17) 

where the slope with respect to r is proportional to 
the average domain size ~ P~^. Hence, the lower the 
temperature, the smaller the kink-antikink creation rate. 
Furthermore, the larger the average metastable size, the 
higher the structural susceptibility to torque. We there- 
fore conclude that the temperature-dependent coopera- 
tive effect is the origin of the observed singularity. 

This is analogous with the instability of a one- 
dimensional (ID) magnetic system. Actually, we see 
the clear correspondence with the ID ferromagnetic Ising 
model: 

T — Tc ^ magnetic field, 
domain wall energy ^ nearest neighbor interaction. 

Hence, it is natural that the role of cooperativity (domain 
wall energy divided by thermal energy) is essential to 
understand the results of Fig. [3l Furthermore, in the 
low temperature limit, there occurs an instability at the 
external torque of r = Tc, which, in the Ising system, 
corresponds to zero external magnetic field |29l |. 

Finally, we shall show that the twist dependence of 
transition rate reported in the experiment of Ref. [3] is 
qualitatively explained by this nucleation picture. Based 
on the above arguments on the nucleation energy (p!3|) , 
the nucleation energy from stable structure and from 
metastable structure are, respectively, given by 

£;nuci = 2(£;b-z + 7k-^c|), 

Erc\^ = 2(^B-Z-7|r-r,|). 

Hence, the transition rates from B to Z structure and 
from Z to B structure, respectively, turn out to be, for 

r > Tc, 

kB-z = t^^e-^^-S kz-B = tz^e-^^-"^ (18) 
and, for r < 

kB-z = ti^e-/'^"- , fcz-B = ^e-''^"- , (19) 



where tB^tz are the inverses of frequency factors for B 
and Z structures. From these results, the equilibrium 
constant Kq^ (= fe-z/^z-e) is given by l^e"^^^^'^"'^^^. 
By noting that r is related to superhelical density a as 
r = Cujqcf i30| and that B-Z transitions were observed at 
(J = (jc ~ —0.01 in the experiment Q, we estimate that 
C = Tc/oooCTc ^ 130 X 10"^°J. Then, we obtain 

Keq = ^e-4/37Cu;o(a-ae)^ (20) 
tB 

where the coefficient of a in the exponent is ^fijCuJo ^ 
4.5 X 10^. Here we used T = 40°C and 7 = 0.37 [fit- 
ted value in Fig. [3fd)]. Thus, the value of the coefficient 
of (J estimated here qualitatively agrees with or, more 
properly speaking, is approximately one-fifth of the ex- 
perimental value (~ 1.4 X 10^) [3l[. 

In summary, we have shown that our nucleation model 
not only provides a clear picture for the mechanically in- 
duced structural transition in DNA as the property of a 
one dimensional chiral material, but also offers an expla- 
nation for recent experimental observations. 



V. SUMMARY 

We have constructed the single DNA mechanical 
model, which describes the interplay between intrinsic 
base-pair structures and global conformations. With this 
model, the mechanical responses of linear DNA to ex- 
ternal torques were simulated and the singular response 
near r = Tc was found in the low-temperature region. 
This singularity gives rise to instability of the B-Z tran- 
sition with minute negative superhelicity, as in the exper- 
iment [3] , which was explained in terms of the cooperative 
effect depending on the average size of the metastable 
domain. Furthermore, the twist dependency of the tran- 
sition rate was qualitatively explained by the nucleation 
picture. These observations are analogous with the in- 
stability of a one-dimensional magnetic system: In low 
temperature limit, there occurs an instability at zero ex- 
ternal magnetic field, which, in our system, corresponds 
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to the external torque of r = Tc [29^. From these esti- 
mates, we have a clear picture for the mechanically in- 
duced structural transition in DNA as one dimensional 
chiral material. 



We expect that this model will provide an additional 
basis for elucidating the structure-configuration interplay 
in higher-order biomolecular architectures, such as nucle- 
osomes. 
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